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We evaluate the observational constraints on the spectral index n, in the context of the ACDM 
hypothesis which represents the simplest viable cosmology. We first take n to be practically scale- 
independent. Ignoring reionization, we find at a nominal 2-a level n ~ 1.0 ± 0.1. If we make the 
more realisitic assumption that reionization occurs when a fraction / ~ 10 -5 to 1 of the matter has 
collapsed, the 2-a lower bound is unchanged while the l-cr bound rises slightly. These constraints are 
compared with the prediction of various inflation models. Then we investigate the two-parameter 
scale-dependent spectral index, predicted by running-mass inflation models, and find that present 
data allow significant scale-dependence of n, which occurs in a physically reasonable regime of 
parameter space. 
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I. INTRODUCTION 

It is generally supposed that structure in the Universe originates from a primordial gaussian curvature perturbation, 
generated by slow-roll inflation. The spectrum Vn(k) of the curvature perturbation is the point of contact between 
observation and models of inflation. It is given in terms of the infiaton potential V(<p) byQ 

where the potential and its derivatives are evaluated at the epoch of horizon exit k = aH. To work out the value of 
(j> at this epoch one uses the relation 

ln(fc cnd /fc) = N(k) = M p 2 / (V/V')dcj>, (2) 



where N(k) is actually the number of e-folds from horizon exit to the end of slow-roll inflation. At the scale explored by 
the COBE measurement of the cosmic microwave background (cmb) anisotropy, -/V(fccoBE) depends on the expansion 
of the Universe after inflation in the manner specified by Eq. ( p0| ) below. 

Given this prediction, the observed large-scale normalization V^ 2 ~ 10~ 5 provides a strong constraint on models 
of inflation. Taking that for granted, we are here interested in the scale-dependence of the spectrum, defined by the, 
in general, scale-dependent spectral index n; 

According to most inflation models, n has negligible variation on cosmological scales so that V-r oc fc n_1 , but we shall 
also discuss an interesting class of models giving a different scale-dependence. 
From Eqs. ([!]) and (§), 

n-l = 2M$(V"/V) - 3M$(V'/V) 2 , (4) 
and in almost all models of inflation, Eq. (|J) is well approximated by 

n - 1 = 2Ml{V"/V) . (5) 



x As usual, M P = 2.4 x 10 18 GeV is the Planck mass, a is the scale factor and H — a/ a is the Hubble parameter, and k/a is 
the wavenumber. We assume the usual slow-roll conditions Mp\V" /V\ -C 1 and Mp(V' /V) 2 <SC 1, leading to 3Htj> ~ —V'. 
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We see that the spectral index measures the shape of the inflaton potential V ((/)), being independent of its overall 
normalization. For this reason, it is a powerful discriminator between models of inflation. 

The observational constraints on the spectral index have been studied by many authors, but a new investigation is 
justified for two reasons. On the observational side, the cosmological parameters are at last being pinned down, as is 
the height of the first peak in the spectrum the cmb anisotropy. No study has yet been given which takes on board 
these observational developments, while at the same time taking on board the crucial influence of the reionization 
epoch on the peak height. On the theory side, it is known that the spectral index may be strongly scale-dependent if 
the inflaton has a gauge coupling, leading to what are called running-mass models. The quite specific, two-parameter 
prediction for the scale dependence of the spectral index in these models has not been compared with presently 
available data. 

II. THE OBSERVATIONAL CONSTRAINTS ON THE PARAMETERS OF THE ACDM MODEL 

Observations of various types indicate that we live in a low density Universe, which is at least approximately flat 



10- 13 1 . In the interest of simplicity we therefore adopt the ACDM model, defined by the requirements that the 
Universe is exactly flat, and that the non-baryonic dark matter is cold with negligible interaction. Essentially exact 
flatness is predicted by inflation, unless one invokes a special kind of model, or special initial conditions. Also, there 
is no clear motivation to modify the cold dark matter hypothesis.^] We shall constrain the parameters of the ACDM 
model, including the spectral index, by performing a least-squares fit to key observational quantities. 

A. The parameters 

The ACDM model is defined by the spectrum Vn(k) of the primordial curvature perturbation, and the four pa- 
rameters that are needed to translate this spectrum into spectra for the matter density perturbation and the cmb 
anisotropy. The four parameters are the Hubble constant h (in units of 100 km s -1 Mpc -1 ), the total matter density 
parameter Oo, the baryon density parameter f2t>, and the reionization redshift zr. As we shall describe, 2r is estimated 
by assuming that reionization occurs when some fixed fraction / of the matter collapses. Within the reasonable range 
/ ~ 10~ 4 to 1, the main results are insensitive to the precise value of /. 

The spectrum is conveniently specified by its value at a scale explored by COBE, and the spectral index n(k). We 
shall consider the usual case of a constant spectral index, and the case of running mass models where n(k) is given 
by a two-parameter expression. Since Vii(kcOBE) is determined very accurately by the COBE data (Eq. ( Jlq ) below) 
we fix its value. Excluding zr and Vn(kcoBE): the ACDM model is specified by four parameters in the case of a 
constant spectral index, or by five parameters in the case of running mass inflation models. 

B. The data 

To compare the ACDM model with observation, we take as our starting point a study performed a few years ago 
H . We consider the same seven observational quantities as in the earlier work, since they still summarize most of the 
relevant data. Of these quantities, three are the cosmological quantities h, f^o, ^b, which we are also taking as free 
parameters. The crucial difference between the present situation and the earlier one is that observation is beginning 
to pin down h and Qq. Judging by the spread of measurements, the systematic error, while still important, is no 
longer completely dominant compared with the random error. At least at some crude level, it therefore makes sense 
to pretend that the errors are all random, and to perform a least squares fit. The adopted values and errors are given 
in Table 1, and summarized below. In common with earlier investigations, we take the errors to be uncorrelated. 

a. Hubble constant On the basis of observations that have nothing to do with large scale structure it seems very 
likely that h is in the range 0.5 to 0.8. We therefore adopt, at notionally the 2-cr level, the value h = 0.65 ± 0.15, 
corresponding to h — 0.65 ± 0.075 at the notional 1-a level. 

b. The matter density The case of the total density parameter Qq is similar to that of the Hubble parameter. 
On the basis of observations that have nothing to do with large scale structure, it seems very likely that Slo lies 
between 0.2 and 0.5, and we adopt at the notional 1-a level the value flo — 0.35 ± 0.075. 



2 In particular, the rotation curves of dwarf galaxies may be compatible with cold dark matter [Q. 
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c. The baryon density As described for instance in [|4|j5j, the baryon density parameter fib has two likely ranges. 
At the l-a level, these are estimated in [Q to be fl h h 2 = .019 ± .002 and fl h h 2 = .007 ± .0015. We adopt the high Q h 
range, which is generally regarded as the most likely, though our conclusions would be much the same if we were to 
adopt the low range. 

d. The rms density perturbation at 8/i _1 Mpc Primarily through the abundance of rich galaxy clusters, a useful 
constraint on the primordial spectrum is provided by the rms density contrast, in a comoving sphere with present 
radius R ~ lOh^ 1 Mpc, at redshift z — to a few. The constrained quantity is conventionally taken to be the present, 
linearly evolved rms density contrast at R = 8/i _1 Mpc, denoted by erg. A recent estimate |J based on low- redshift 
clusters gives at l-a 

a 8 = a 8 ^ QA7 (6) 
fjg = .560 ± .059 . (7) 

This constrains the primordial curvature perturbation on the scale k ~ k s = Mpc) -1 . 

e. The shape parameter The slope of the galaxy correlation function on scales of order lh^ 1 to 100/i _1 Mpc is 
conveniently specified by a shape parameter || T, defined by 

f = r-0.28(n 8 1 -l) (8) 
F = Q hexp(-Q. B -Q B /ilo)- (9) 

(The quantity T determines, to an excellent approximation, the shape of the matter transfer function on scales /c _1 ~ 1 
to lOO/i -1 Mpc, while the second term accounts for the scale dependence of the primordial spectrum. For definiteness, 
we evaluate n at k — ks, in the case that n has significant scale dependence.) A fit reported in [0] gives T = .23 with a 
15% uncertainty at 2-er. A more recent fit with more data JlJ gives T — .20 to .25, depending on the assumed velocity 
dispersion, but with 15% statistical uncertainty at the l-a level.| We therefore adopt V = .23, with 15% uncertainty 
at l-a. 

f. The COBE normalization of the spectrum To a good approximation, the spectrum Ct of the cmb anisotropy at 
large £ is sensitive to the primordial spectrum on the corresponding scale at the particle horizon, 

k(£,n ) = — e — (io) 

x hol = 2H- 1 n~ 1/2 (1 + 0.084 In Q ) . (11) 

The COBE measurements cover the range 2 < £ < 30, and they constrain Vnik) on the corresponding scales. Instead 
of Viz, it is usual in this context to consider a quantity 5h , which is of direct interest for studies of structure formation 
and is defined by 

Mfc)4^W /2 w (12) 

5 S io 

gN ,^fl + ^-l + ^)". (13) 

The factor g/£lo, normalized to 1 at f2o = 1, represents the f2o - dependence of the present, linearly evolved, density 
contrast after pulling out the scale-dependent transfer function and Vtz- Equivalently, a(Q)g(fl) is the time-dependence 
of the density contrast after matter domination. 

According to the ordinary (as opposed to 'integrated') Sachs- Wolfe approximation 

Att r°° (\k 

Ce =25j -y j2e{kXhol)Vnik) - (14) 

In the regime £ 1, it satisfies Eq. ( |ic| ) because jf peaks when its argument is equal to £. In the ACDM model, 
the Sachs- Wolfe approximation is quite good in COBE regime, but still the quality of the data justify using the full 
(linear) calculation, given for instance by the output of the CMBfast package |Tq|. 



3 See Table 3 of [H; in the present context one should focus on the last three rows of the Table. 
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Consider first the case n = 1 (scale- independent spectrum). In the Sachs- Wolfe approximation, the value of V-r 
obtained by fitting the COBE data is independent of the cosmological parameters h, Slo and ^t>- Using instead the 
full calculation, a fit to the data by Bunn and White || gives 

X o-0-785-0.051nn T 

oh = "o "h 
10 5 fe = 1.94 ±0.08, (15) 



As expected, the corresponding spectrum of the curvature perturbation has only mild dependence on Qq CPn oc 

Consider next the case of a scale- independent spectral index n =/= 1. Dropping an insignificant term quadratic in 
n — 1, the fit of Bunn and White B handles the n-dependence by assuming that Eq. (|l^) holds at a 'pivot' scale 
&cobe which is independent of f2o □ 

&cobe = 6.6-ffo , ■ (16) 

Insofar as the approximation Eq. ( |l0| ) is valid, this corresponds to fixing Ce at an r2o-dependent value of £, which is 
£ = 13 for n , and I = 22 for our central value Sl = -35. 

In the case of a scale- independent n, an alternative fit is provided by the CMBfast package, which chooses VvSk) 
to fit an n-independent best-fit value of Ciq. As expected, the output of CMBfast is in good agreement with the 
Bunn- White fit. Even better agreement is obtained using 



«COBE 



(fi ) = 13.2/x hor , (17) 



which reduces to Eq. ([IB]) for Qq = 1. Insofar as Eq. ( |l0| ) is valid, this fio-dependent pivot for k corresponds to an 
f^o-independent pivot for £, namely £ = 13. 

We are also interested in the scale-dependent n predicted by the running-mass inflation models. However, as the 
range of scales explored by COBE corresponds to only AN ~ 2, with the central values of £ the most important, wc 
can take the variation of n to be negligible on these scales. 

Guided by these considerations, we have adopted three slightly different versions of the COBE normalization, chosen 
for convenience according to the context. When calculating T and a$, we in all cases fixed 6h at the central value 
given by Eq. ([l5]), at the Bunn- White pivot point fccoBE- When calculating the height of the first peak in the cmb 
anisotropy, in the case of the running-mass model, we used Eq. (|To|), with Sh again fixed at the central value given by 
Eq. ( |l5| ) but now evaluated at the slightly more accurate pivot point /ccobe(^o)- Finally, when evaluating the peak 
height in the case of scale-independent n, we used a linear fit to the output of CMBfast. Explicit expressions for the 
peak height will be given after considering the effect of reionization. 

g. The peak height The model under consideration predicts a peak in the cmb anisotropy at £ ~ 210 to 230, 
and presently available data JlOj— jT3[| confirm the existence of a peak at about this position. We adopt as a crucial 
observational quantity C pea k, defined as the maximum value of 

C e = £(£+l)Ci/2Tr. (18) 



Presently available data give conflicting estimates [|10|-[13[ of y C pG ak, with central values in the range 70 to 90 /iK. 
We adopt (80 ± 10) fj,K with the uncertainty taken to be at 1-a. 



C. Reionization 



The effect of reionization on the cmb anisotropy is determined by the optical depth r. We assume sudden, complete 
reionization at redshift zr,, so that the optical depth r is given by |lj,|l5| 



r = 0.035^ 

"0 



/i(V^o(l + 2 R ) 3 + l-r!o-l) • (19) 



4 Keeping the quadratic term, the 'pivot' scale at which Eq. (^5|) holds is dependent on n, but still independent of flo. A 
related fit by Bunn, Liddle and White J9| keeps a cross-term in (n — 1) and fio, which makes the 'pivot' scale increase with 
CIq 1 , though not as strongly as in Eq. (hj) below. 
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In previous investigations, zr has been regarded as a free parameter, usually fixed at zero or some other value. 
In this investigation, we instead take on board that fact that zr can be estimated, in terms of the parameters that 
we are varying plus assumed astrophysics. Indeed, it is usually supposed that reionization occurs at an early epoch, 
when some fraction / of the matter has collapsed, into objects with mass very roughly M — 10 6 M©. Estimates of / 
are in the range |i6| ] 

l(T 4 - 4 </<l. (20) 
In the case /< 1, the Press-Schechter approximation gives the estimate 

l + ^^M^erfc-V) (/«!)• (21) 

Here a(M) is the present, linearly evolved, rms density contrast with top-hat smoothing, and 5 C = 1.7 is the overdensity 
required for gravitational collapse. ( g is the suppression factor of the linearly evolved density contrast at the present 
epoch, which does not apply at the epoch of reionization.) In the case / ~ 1, one can justify only the rough estimate 

1 + (/-I)- (22) 

(This estimate is not very different from the one that would be obtained by using / = 1 in Eq. (pi|).) 

In our fits, we fix / at different values in the above range, and find that the most important results are not very 
sensitive to / even though the corresponding values of zr can be quite high. 



D. The predicted peak height 

The CMBfast package jL8| gives Cg, for given values of the parameters with n taken to be scale-independent. 
Following , we parameterize the CMBfast output at the first peak in the form 

220 y /2 

10" ) 



Cpcak - \j Cpelk ( TFT ) > (23) 



where 

v = a n {n - 1) + a h ]xi(h/0M) + a ln(fi /0.35) + a h h 2 (Q h - f2 b (0) ) - 0.65/ (t)t . (24) 



Cpeik i s the value of y C pG ak evaluated with each term of v equal to zero. The coefficients for the high choice 

n { ° ] h 2 = 0.019 are a n = 0.88, a h = -0.37, a = -0.16, a h = 5.4, and = 77.5 ^K. The formula reproduces the 

CMBfast results within 10% for a 1-er variation of the cosmological parameters, h, flo and fib, and ncoBE = 1.0±0.05. 



With the function /(r) set equal to 1, the term — 0.65t is equivalent to multiplying y C poa k by the usual factor exp(— r). 
We use the following formula, which was obtained by fitting the output of CMBfast, and is accurate to a few percent 
over the interesting range of r; 

/ = l-0.165r/(0.4 + r). (25) 

For the running-mass model, we start with the above estimate for n = 1, and adjust it using Eq. (|To|). Adopting 
the COBE normalization mentioned earlier, this adjustment is 

_ 5 H (k(£,n )) 



V °poak 



In the case of constant n, this prescription corresponds to the previous one with a n — 0.91, in good agreement with 
the output of CMBfast. 
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FIG. 1. The nominal 1- and 2-a bounds on n. In the left-hand panel, the reionization redshift is fixed. In the right-hand 
panel, reionization is assumed to occur when a fixed fraction / of matter collapses (corresponding reionization redshift, not 
shown, is roughly in the range 10 to 35). A result n > 1 would rule out most known models of inflation, a result n > .93 would 
rule out 'new' inflation with a cubic potential; these cases are indicated by horizontal lines. 

III. CONSTANT SPECTRAL INDEX 

A. The observational constraints 

Most models of inflation make n roughly scale- independent, over the cosmologically interesting range. We therefore 
begin by considering the case that n is exactly scale-independent. The resulting bound on n is shown in Figure 1. 
In the left-hand panel we make the traditional assumption that reionization occurs at some fixed redshift zr. In the 
right-hand panel we make the more reasonable assumption, that it occurs when some fixed fraction / of the matter 
collapses, in a reasonable range 10 -4 ' 5 < / < 1. The bounds in the latter case are relatively insensitive to /, because 
the corresponding range of zr is narrower; everywhere on the displayed curves, zr is within (usually well within) the 
range 8 to 36. Details of the fit for zr = 20 are given in Table 1. Practically the same fit is obtained if instead we fix 

/ at ltr 1 - 9 . 

The least-squares fits were performed with the CERN minuit package, and the quoted error bars invokes the usual 
parabolic approximation (i.e., it they are the diagonal elements of the error matrix). The exact error bars given by 
the same package agree to better than 10%. For zr, our results are similar to those obtained in [^o|, but more precise 
because of improvements in our knowledge of the cosmological parameters; they are also similar to those obtained in 
pjj , if we take the errors to be the ones given by the error matrix. (We do not know why the exact error bars in pi]] 
are about three times bigger, in conflict with both our work and that of P0|.) 

After we completed this work, the BOOMeranG JnJ and MAXIMA p3| measurements of the cmb anisotropy 
appeared, both of which extend to the second acoustic peak. Fits to these data (3J]j3(| seem to again give a similar 
constraint on n, but the values for 17b, and h outside our adopted 2-a range. At the time of writing, the new cmb 
data have not been included in a fit of the type that we are performing (i.e., with with strong prior requirements on 
the cosmological parameters, as well as on the small-scale data as and T). 

B. Models of inflation giving n < 1 

Although the quality and quantity of data are insufficient for a proper statistical analysis, these bounds on n are 
very striking when compared with theoretical expectations. These expectations ]22| , |l"5| ] are summarize of] in Tables 



5 This Table excludes the running-mass models to be discussed later, and a recently-proposed model |i7| giving (n — l)/2 = 
— 2/N. It also excludes the ad hoc 'chaotic inflation' potentials V oc (f> p , which give n — 1 = —(2 +p)/(2N) with a significant 
gravitational contribution to the cmb anisotropy. 
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0.35 


0.65 


0.23 
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80 fjK 


error 




0.002 


0.075 


0.075 


0.035 


0.059 


10 imK 


fit 


1.064 


0.019 


0.34 


0.63 


0.19 


0.59 


77 /iK 


error 


0.077 


0.002 


0.06 


0.06 








x 2 




9 x 10 -5 


3 x 10~ 2 


0.1 


1.3 


0.2 


0.1 



TABLE I. Fit of the ACDM model to presently available data, with zr = 20. The spectral index n is a parameter of the 
model, and so are the next three quantities. Every quantity except n is a data point, with the value and uncertainty listed in 
the first two rows. The result of the least-squares fit is given in the lines three to five. All uncertainties are at the nominal 1-cr 
level. The total \ 2 is 1-8 f° r t wo degrees of freedom. 



TABLE II. Predictions for the spectral index n(k). Wavenumber k is related to number of e-folds N by dlnfc = —dN. 
Constants q and Q are positive, and p can have either sign. 



Comments 


v{4>)/v 


i(n- 


1) 


Mass term 




±M$m 2 /V 


p integer < — 1 or > 3 


1 + |c|0 p 


p-i 

p-2 N 


1 

max-iV 


Spont. broken susy 


l + |c|In£ 


1 

2N 




Various models 


1 - e -9 * 


1 

N 




p > 2 or — oo < p < 1 


1 - \c\(t> p 


-(*= 

VP- 


2} N 



2 and 3, and we now discuss them beginning with the usual case n < 1 (red spectrum), 
references are given in |p2|| . 

The simplest prediction is for a potential of the form[] 

V = Vn- \m 2 <b 2 + • • • , 



Details of the models and 



(27) 



leading to n — 1 = — 2Mpm 2 /Vq. This is the form that one expects if is a string modulus (Modular Inflation), or a 
pseudo-Goldstone boson (Natural Inflation) , or the radial part of a massive field spontaneously breaking a symmetry 
(Topological Inflation). The vacuum expectation value (vev) of <fi in these models is expected to be of order Mp or less, 
while the potential Eq. ( |27| ) gives ((f)) ~ (1 — n)~ x l 2 Mp. Therefore, the present bound n > 0.9 is already beginning 
to disfavor these models. The potential Eq. ( p7| ) may however give n very close to 1 if the potential steepens after 
cosmological scales leave the horizon, for instance in an inverted hybrid inflation model. 

Of the remaining models of Table 2, those giving a red spectrum involve a potential basically of the form 



V = V (1 + c(j? + ■ ■ ■) 



(28) 



with c negative and p not in the range 1 < p < 2. ( 'New' inflation corresponds to p an integer > 3, while mutated 
hybrid inflation models account for the rest of the range. The logarithmic and exponential potentials in Table 2 may 
be regarded as the limits respectively p — > and p — > — oo.) With this form, the prediction is 



1 



p-1 
p-2 



2 

N 



(29) 



6 In this expression and in Eqs. (|28|) and (]35|), the remaining terms are supposed to be negligible, and Vb is supposed to 
dominate, while cosmological scales leave the horizon. 
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TABLE III. Predictions for the spectral index n, for some potentials of the form Vb(l + C(f> p ) with negative c. The case p — > 
corresponds to the potential Vb(l + cln^), and the case p — » — oo corresponds to Vb(l — e~ 91 *). The parameter -/Vcobe < 60 
depends on the cosmology after inflation. 



p 


/i 






AfcoBE = 50 


iVcOBE = 20 


p -> 


0.98 


0.95 


p = -2 


0.97 


0.93 


p — ► ±oo 


0.96 


0.90 


p = 4 


0.94 


0.85 


p = 3 


0.92 


0.80 



For the moment, we ignore the mild scale-dependence and set TV = Acobe- 
Depending on the history of the Universe, 

A^cobe * 60 - ln(10 16 GeV/U 1 / 4 ) - i ln^ 1 / 4 /T reh ) - N . (30) 

In this expression, T rc h is the reheat temperature, while the final contribution —No (negative in all reasonable cos- 
mologies) encodes our ignorance about what happens between the end of inflation and nucleosynthesis. Let us pause 
to discuss this ignorance. In the present context, we are defining T ro h as the temperature when the Universe first 
becomes radiation dominated after inflation. In the conventional cosmology, radiation domination persists until the 
present matter dominated era begins, long after nucleosynthesis. If this is the case, and if also slow-roll inflation 
gives way promptly to matter domination as is the case in most models, then iVo = 0.Q In this conventional case, 
Acobe is largely determined by Uq 1 ^ 4 , and hence by the model of inflation. It is certainly in the range 32 to 60 (lower 
limit corresponding to Vq^ 4 = 100 GeV) and much more likely in the range 40 to 60 (lower limit corresponding to 

U 1/4 ~ 10 10 GeV and T rch ~ 100 GeV). 

However, the conventional cosmology need not be correct. In particular, the initial radiation-dominated era may 
give way to matter domination by a late-decaying particle, and most crucially there may be an era of thermal inflation 
p3j during the transition. This unconventional cosmology, with its huge entropy dilution after inflation, is indeed 
demanded in many inflation models, if gravitinos created from the vacuum fluctuation [p4f persists to late times [p5| . 
Even one bout of thermal inflation will give iVo ~ 10 and additional bout(s) cannot be ruled out. Thus, from the 
theoretical viewpoint, -/Vcobe can be anywhere in the range to 60. 

Let us discuss the prediction Eq. ( p9] ) , excluding for simplicity the ranges < p < 1 and 2 < p < 3 (recall that the 
straightforward 'new' inflation models make p an integer > 3). Taking the maximum value -/Vcobe — 60, we learn 
that n < 0.93 for p = 3 (the lowest prediction), and n < 0.95 for p = 4. Looking at the right-hand panel of Figure 
1, we see that at nominal 1-a level, the former case is ruled out, though it is still allowed at the 2-a level. Stronger 
results hold in the if -/Vcobe < 60. Looking at things another way, a lower bound on n gives a lower bound on A^cobe, 

Ncobe > V —^z ■ (31) 

p — L 1 — n 

Even with present data, the 2-cr result n > .9 gives A^cobe > 40 for p = 3, and A/cobe > 20 for p> 3. 
The scale dependence given by Eq. (p9^is 

(n-l) 2 <0. (32) 



dn 1 / p 



dlnfc 2 \p- 1 

Over the cosmological range of scales m(£;/£;coBE) ^ s a * most a few, and in particular ln(8 _1 /iMpc _1 //ccoBE) — 4, 
corresponding to 



7 In some inflation models, slow-roll is followed by an extended era of fast-roll giving No of order a few; for simplicity we ignore 
that possibility in the present discussion. 
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X X2 



An = n s - ncoBE = -.02 I ^— -J [~qY J <°- ( 33 ) 

Taking n = 0.9 to saturate the present bound, this gives |An| < 0.02 with p > 3, and |An| < 0.04 with p < 0. Even 
in the latter case, the change in n is hardly significant with present data. 

C. Models giving n > 1 

Known models giving n > 1 (blue spectrum) are all of the hybrid inflation type. The simplest case is V = 
Vo + \vr?(\) 2 \ it gives the scale-independent prediction n — 1 = 2A/pm 2 /Vb, which may be either close to 1 or well 
above 1. 

The other cases involve a potential of the form V = Vo (1 + c(j) p ) with positive c, and p an integer > 3 or < — 1. 
There is a maximum (early-time) value for N, and the prediction 

n-l = P —l l - . (34) 

Barring the fine-tuning iVcoBE — JV max , this gives n — 1 <C 0.04, which is compatible with the observational bound. 
The scale-dependence of n in these models is still given by Eqs. ([32|) and ([53|); it may be observationally significant 
only in the fine-tuned case Nqobe — JV max , which we have not investigated. 

IV. THE RUNNING MASS MODELS 
A. The potential 

We have also done fits with the scale-dependent spectral index, predicted in inflation models with a running inflaton 
mass p7|-|32|. In these models, based on softly broken supersymmetry, one-loop corrections to the tree-level potential 
are taken into account, by evaluating the inflaton mass-squared m 2 (ln(Qj) at the renormalization scale Q ~ 0,^ 

V = Vv + ^m 2 (\n{Q))<f + ■■■ . (35) 

Over any small range of 0, it is a good approximation to take the running mass to be a linear function of In </>. 
This is equivalent to choosing the renormalization scale to be within the range, and then adding the loop correction 
explicitly, 



1 ,„ _ 1 „ ^Vo_ 



V = V + -m 2 {hxQ)ct> 2 - ~ c (lnQ)^ 2 ln^/Q) . (36) 



The dimensionless quantity c specifies the strength of the coupling. Let us discuss its likely magnitude, taking for 
definiteness Q — 0cobe- 

It has been shown fl3Cfl that the linear approximation is very good over the range of <f> corresponding to horizon exit 
for scales between &cobe anc ^ 8ft. -1 Mpc. We shall want to estimate the reionization epoch, which involves a scale of 
order k~ eion ~ 10 -2 Mpc (enclosing the relevant mass of order 10 6 M Q ). Since only a crude estimate of the reionization 
epoch is needed, we shall assume that the linear approximation is adequate down to this 'reionization scale'. In other 
words, we assume that it is adequate for between </>cobe and re ion, the subscripts denoting the value of 4> when 
the relevant scale leaves the horizon. Within this range, we it is convenient to write Eq. ([57]) in the form J30[ ] 

1 Vo_ ,o (, 4> l 



8 The choice Q ~ <j> is to be made in the regime where (f> is bigger than the relevant masses. When Q falls below the relevant 
masses, m 2 (Q) becomes practically scale-independent (the mass 'stops running'). We have a running mass model if inflation 
takes place in the former regime, which happens in some interesting cases , including that of a gauge coupling. 
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so that 



In these expressions, the constants c and <f>* both depend on the renormalization scale Q, which can be chosen anywhere 
in the range corresponding to cosmological scales (say Q = </>cobe)- The dimensionful constant is related to the 
mass-squared by 

ln(0*/Q) = — — — -. (39) 

c(Q)V 2 

Note that the limit of no running, c — > 0, corresponds to finite c| \sx{4> / (f>*)\, so that Eq. ( p7[ ) in that limit gives back 
Eq. ( p5| ) with a constant mass. 

In general, the point <p = cj>* may be far outside the regime where the linear approximation Eq. (|37| ) applies. 
However, in simple models the cosmological regime is sufficiently close to that point that the linear approximation is 
approximately valid there. In that case, we can trust the Eq. ( |37|) and its derivatives for <f> = 0*; since V' vanishes at 
that point, there are four clearly distinct models of inflation as shown in Figure §. The labeling (i), (ii), (iii) and (iv) 
is the one introduced in [po| . In Models (i) and (ii), c is positive and the potential has a maximum near 0„, while in 
Models (iii) and (iv), c is negative and there is a minimum. In Models (i) and (iv), <fi moves towards the origin, while 
in Models (ii) and (iii) the opposite is true. Even if Eq. ( |37| ) is not valid near <j> = <fi*, this fourfold classification of 
models, according to the sign of c and the direction of motion of (f>, is still useful. 

Let us discuss the likely magnitude of c, assuming that a single coupling dominates the loop correction. The value 
of c is conveniently obtained from the well-known RGE for dm 2 /d(lnQ). If a gauge coupling dominates one finds Eq] 



V c 2C . , 



am 1 . (40) 

p " 



Here, C is a positive group-theoretic number of order 1, a is the gauge coupling, and m is the gaugino mass. We see 
that if the loop correction comes from a single gauge coupling, c is positive, corresponding to Model (i) or Model (ii) . 
If a Yukawa coupling dominates, one finds plj| (for negligible supersymmetry breaking trilinear coupling) 



Vac D , 2 2 



wr-^ lxlmLp ' (41) 

where D is a positive constant counting the number of scalar particles interacting with the inflaton, ?7ij 2 oop is their 
common susy breaking mass-squared, and A is their common Yukawa coupling. In this case, c can be of either sign. 

To complete our estimate of c, we need the gaugino or scalar mass. The traditional hypothesis is that soft super- 
symmetry breaking is gravity-mediated, and in the context of inflation this means that the scale Ms of supersymmetry 
breaking will be roughly V^ 4 . (As usual we are defining Ms = y/F, where F is the auxiliary field responsible for 
spontaneous supersymmetry breaking in the hidden sector. We also assume that there is no accurate cancelation in 
the formula V = \F\ 2 — 3Mpiny 2 , which is the case in most supersymmetric inflation models p2|.) With gravity- 
mediated susy breaking, typical values of the masses are m 2 ~ |-w 2 | ~ Vo/Mp, which makes |cj of order of the 
coupling strength a or |A| 2 . At least in the case of a gauge coupling, one then expects 

|c| - 10" 1 to 10~ 2 . (42) 

In special versions of gravity-mediated susy breaking, the masses could be much smaller, leading to |c| <C 1. In that 
case, the mass would hardly run, and the spectral index would be practically scale-independent. With gauge-mediated 
susy breaking, the masses could be much bigger; this would not lead to a model of inflation (unless the coupling is 
suppressed) because it would not satisfy the slow-roll requirement |c| < 1. 



B. The spectrum and the spectral index 

Using Eq. (p|) we find 



se cAN(k) = cln (^) (43) 

AA^(fc) = Ncobe - N(k) = ln(fc/fc C OB E ) , (44) 
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Model (i) 



Model (ii) 




Model (iii) 




FIG. 2. Sketches of the potential for the different models in the case an extremum exists: the right panel shows the inflaton 
behavior for Models (i) and (ii), while the left panel shows Models (iii) and (iv). 



where s is an integration constant |] Eq. (g) then gives 



i(k) - 1 



= se cAN ^-c. 



(45) 



Some lines of fixed ncoBE in the plane s versus c are shown in the left-hand panel of Figure [| In order to evaluate 
Eq. (Eq), we also need the variation of Sh which comes from integrating this expression, 



S H (k) 



fo(fccOBE) 



exp 



„cAN 



1 



cAN 



(46) 



We are mostly interested in cosmological scales between /ccobe and kg, corresponding to < A A" < 4. In this 
range the scale-dependence of n is approximately linear (taking |c| < 1) and the variation An = n§ — ticobe is given 
approximately by 



dlnfc 



(47) 



In contrast with the prediction Eqs. ( p2|) and (|33|) of the earlier models we considered, An is positive. Also in contrast 
with those models, it is not tied to the magnitude of |n— 1|, and (as we shall see) may be significant even with present 
data, for physically reasonable values of the parameters. In the right-hand panel of Figure [|, we show the branches 
of the hyperbola 8sc = An, for the reference value An = 0.04. Within the hyperbola, the scale-dependence of n is 
probably too small to be significant with present data. 

The spectral index Eq. ( [45| ) depends on the coupling c, which we already discussed, and the integration constant 
s. To satisfy the slow-roll conditions Mp\V" /V\ <C 1 and Mp(V /V) 2 <C 1, both c and s must be at most of order 
1 in magnitude. Significant additional constraints on s follow, if we make the reasonable assumptions that the mass 
continues to run to the end of slow-roll inflation, and that the linear approximation remains roughly valid. Indeed, 
setting /S.N — A^cobe, Eq. (^|) becomes s — e~ cArcOBE cln(0»/0 cn d). Discounting the possibility that the end of 
inflation is very fine-tuned, to occur close to the maximum or minimum of the potential, this gives a lower bound 



\s\ > e - cNcoBE 



c . 



(48) 



In the case of positive c (Models (i) and (ii)), we also obtain a significant upper bound by setting AA^ = Aqobe in 
Eq. (f45|), and remembering that slow- roll requires |n — 1| < 1; 



s <e 



-cJVcOBE 



(c>l). 



(49) 



In an earlier paper pcj we used a = se 



cJVcOBE 



but s is more convenient. 
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In the simplest case, that slow-roll inflation ends when n — 1 actually becomes of order 1, this bound becomes an 
actual estimate, \s\ ~ 6 - cN cobe_ 

In the case of Models (i) and (iv), the mass may cease to run before the end of slow-roll inflation (but after 
cosmological scales leave the horizon, or the running mass model would not apply) at some point N Tnn . In this 
somewhat fine-tuned situation, Acobe in the above estimates should be replaced -/Vcobe — A run , which may be much 
less than -ZVcobe- In the case of Model (iv), this leads to a weaker lower bound 

s > |c| (c < 0) . (50) 

In the case of Model (i) it leads to a weaker upper bound 

s < 1 (c> 0) . (51) 

In the left-hand panel of Figure a, we show the bounds relevant to the choice of parameter's ranges, i.e. the lower 
bound Eq. (EJj) , the upper bound Eq. (EJJ) and the weak lower bound Eq. (EG) . 



C. The magnitude of the spectrum 

Although it is not directly relevant for our investigation of the spectral index, we should mention the constraint on 
the running mass model that comes from the observed magnitude V-r 1 ^ 2 ~ 10~ 5 of the spectrum. From Eq. (Q), 

4^ V Q (2s \ 1 

25 61NL 



^ = 72Xr2 ex P - )U2- ( 52 ) 



This prediction involves V$ and </>*, in addition to the parameters c and s that determine the spectral index. 

The simplest thing is to again assume gravity-mediated susy breaking, with the ultra-violet cutoff at the traditional 
scale around Mp, and the same supersymmetry breaking scale during inflation as in the true vacuum so that Vq^ 4 ~ 
10 10 GeV. In this scenario, one expects |m 2 (Q)| ~ Vo/Mp at Q ~ Mp. As Stewart pointed out in the first paper on 
the subject, with this very traditional set of assumptions, Eq. (|5^) can give the correct COBE normalization, with |c| 
in the physically favored range 10" 1 to 1O~ 2 .0 

It is remarkable that the most traditional set of assumptions can give a model with the correct COBE normalization, 
and, as we shall see, with a viable spectral index. If one relaxes these assumptions, there is much more freedom in 
choosing Vq and (/>*. Such freedom may be very welcome, in coping with the difficulty of implementing inflation in 
the context of large extra dimensions p3| . 

D. Observational constraints on the running mass models 

Extremizing with respect to all other parameters, we have calculated \ 2 m the s vs. c plane and obtained contour 
levels for \ 2 equal to the minimum value plus 2.41 and 5.99 respectively, corresponding nominally to the 70% and 
95% confidence level in two variables. (The x 2 function presents actually two nearly degenerate minima in the allowed 
region, one in the positive and one in the negative quadrants (Models (i) and (iii)), separated by a very low barrier, 
but we assume that the usual quadratic estimate of the probability content is not very far from the true value.) 

The allowed region is shown in the right-hand panel of Figure pL for the case that reionization occurs when / ~ 1. 
For c = or s — the constant n result is recovered with n — 1 — —2c or 2s; our plots give in this case a slightly 
larger allowed interval with respect to the two sigma value in the previous section, due to the mismatch between the 
statistical one variable and two variables 95% CL contours. This allowed region is not too different from the one that 
we estimated earlier by imposing the crude requirement |n— 1| < 0.2 at both the COBE scale and the low scale 
corresponding to Ncobe — 10. (Note that in the earlier work we used the less convenient variable a = s exp(cAcoBE), 
instead of s.) 

The allowed region for Models (ii) and (iv) lies inside the hyperbola corresponding to An = .04, which means that 
their scale-dependence is hardly significant at the level of present data. In contrast, the allowed region for Models 



10 At the crudest level, one can verify this using the linear approximation Eq. ( ^t| ) all the way up to the <p ~ Mp, corresponding 
to ln(Mp/0*) ~ 1/c ~ 10 to 100. Proper calculations p3-K0J using the RGE's lead to the same conclusion. 
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FIG. 3. The parameter space for the running mass model. In the left-hand panel we show the straight lines corresponding 
to ncoBE = 1.2, 1.0 and 0.8. Also shown in the left-hand panel are the lower bound Eq. (^), the upper bound Eq. ([l9|), and 
(diagonal line in upper right quadrant) the weak lower bound Eq. (jib]). (The weak upper bound Eq. (pl[ ) is off the scale.) As 
explained in the text, these curves define the theoretically reasonable region of the parameter space. In the right-hand panel, 
we show the region allowed by observation, in the case that reionization occurs when / ~ 1. Note that the allowed region is 
parallel to the fixed tlcobe lines around tlcobe — 1, as one would expect. To show the scale-dependence of the prediction for 
n, we also show in this panel the branches of the hyperbola 8sc = An = ng — kcobe, for the reference value An = 0.04. 



(i) and (iii) extends to An > 0.2, representing an extremely significant scale-dependence even with present data. To 
demonstrate this, we show in Figures |^ and || the allowed regions for Models (i) and (iii) in the n% versus ticobe 
plane. In the case of Model (iii), the theoretical bounds on the parameters restrict the parameter space to a small 
corner of the allowed region, within which n has negligible variation. In contrast, there is no significant theoretical 
restriction on the parameters in the case of Model (i), and n has significant variation in a physically reasonable regime 
of parameter space. In both cases, a lower value of the fraction of collapsed matter / just reduces the allowed region 
at large n, without affecting significantly the allowed scale-dependence of n. 

In the case of Model (i) , a further observational constraint comes from the requirement that the density perturbation 
on scales leaving the horizon at the end of inflation, should be small enough to avoid dangerous black hole formation. 
The linear approximation is not adequate on such small scales, and one should instead evaluate the running mass 
using the RGE. The simplest assumption is that the RGE corresponds to a single gauge coupling, either with or 
without asymptotic freedom The black hole constraint has been evaluated for these cases jj4|. The constraint 
amounts more or less to an upper bound on ncoBE, typically in the range 1.1 to 1.3 depending on the choices of 
Nqobe and other parameters. Such a bound significantly reduces the allowed region of parameter space, but still 
leaves a region where n has a strong variation. 



V. CONCLUSION 



In the context of the ACDM model, we have evaluated the observational constraint on the spectral index n(k). This 
constraint comes from a range of data, including the height of the first peak in the cmb anisotropy, which we take to 
be 80 ± 10/iK (nominal l-er). Reionization is assumed to occur when some fixed fraction / of the matter collapses, 
and the most important results are insensitive to this fraction in the reasonable range 10 -4 < / < 1. 

We first considered the case that n has negligible scale dependence, comparing the observational bound with the 
prediction of various models of inflation. A significant improvement in the 2-er lower bound, which may well occur 
with the advent of slightly better measurements of the cmb anisotropy, will become a serious discriminator between 
models of inflation. Even the present bound has serious implications if, as is very possible, late-time gravitino creation 
or some other phenomenon requires an era of thermal inflation after the usual inflation. 

We also considered the running mass models of inflation, where the spectral index can have significant scale- 
dependence. Because of this scale dependence, it is in this case crucial to fix not the epoch of reionization, but the 
fraction / of matter that has collapsed at that epoch. We presented results for the choice / = 1 (corresponding to 
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FIG. 4. Allowed region in the iicobe — 1 vs ng — 1 plane at 95% CL (solid line) and 70% CL (dashed line) for positive s 
and c (Model (i)). The two panels correspond to different hypotheses about the reionization epoch. In the right panel, it is 
assumed that reionization occurs when a fraction / = 10~ 2,2 of the matter has collapsed into bound structures, while in the 
left panel the fraction is taken to be / ~ 1. 



zr ~ 13 if the spectral index has negligible scale-dependence), and for a perhaps more reasonable choice / = 10 . 
In the running-mass models, the scale-dependent spectral index n(k) is given by n — 1 = sexp(cAiV) — c, where 
AN = ln(fccoBE/fc)- The parameters in this expression can be of either sign, leading to four different models of 
inflation. Barring fine-tuning, one expects s to be in the range |c|e _cArcOBE < \s\ < e _cArcOBE . The parameter c 
depends on the nature of the soft supersymmetry breaking, but in the simplest case of gravity-mediation it becomes 
a dimensionless coupling strength, presumably of order 10 _1 to 10 -2 in magnitude. 

Without worrying about the origin of the parameters c and s, we have investigated the observational constraints 
on them. In the case c, s > (referred to as Model (i)) we find that n can have a significant variation on cosmological 
scales, with n— 1 passing through zero signaling a minimum of the spectrum of the primordial curvature perturbation. 
In a future paper, we shall exhibit the possible effect of this scale-dependence on the cmb anisotropy, at and above 
the first peak. 
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